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Recent literature on delocalization in non-Hermitian systems has stressed criteria based on sensi- 
tivity of eigenvalues to boundary conditions and the existence of a non-zero current. We emphasize 
here that delocaUzation also shows up clearly in eigenfunctions, provided one studies the product of 
left- and right-eigenfunctions, as required on physical grounds, and not simply the squared modulii 
of the eigenfunctions themselves. We also discuss the right- and left-eigenfunctions of the ground 
state in the delocalized regime and suggest that the behavior of these functions, when considered 
separately, may be viewed as "intermediate" between localized and delocahzed. 
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I. INTRODUCTION 

A delocalization phenomenon in a particu- 
larly simple class of non-Hermitian random sys- 
tem has attracted considerable attention re- 
cently @,|,||||,g,|,|jio|jl|ji,0,0jll,§,0§|jl| . 
Among the more recent work is a report by Silve- 
strov [pO|, based on an analysis of eigenfunctions, which 
claims that the phenomenon studied was not actually 
delocalization, but "localization of a very unusual kind." 
Although Silvestrov subsequently revised his views, he 
still maintains that "the transition from real to complex 
spectra in Id disordered systems with (an) imaginary 
vector potential is not a delocalization transition" pl| . 
In this paper, we review some basic facts of the non- 
Hermitian delocalization (Sec. and then take issue 
with Silvestrov's interpretation. We stress in Sec. _|ll| 
that the criteria for delocalization used in Refs. [^P] 
are entirely consistent with a conventional one based on 
eigenfunctions, provided one studies the correct physical 
quantity, namely the product of the left- and right- 
eigenfunctions associated with a given state. In Sec. EV, 
we comment on the interesting results of Silvestrov for 
left- and right-eigenfunctions considered separately for 
large asymmetry parameter. We show for the ground 
state that the results are related to earlier results ob- 
tained for charge density waves |^ and population biol- 
ogy [T^ . From this viewpoint, we argue that the behavior 
of the left- and the right-eigenfunctions is "intermediate" 
between localized and delocalized behavior. 



II. NON-HERMITIAN DELOCALIZATION: 
EIGENVALUES AND CURRENT 

Let us first review some basic facts about non- 
Hermitian delocalization. A typical example of the sys- 
tems in question is the one-particle Hamiltonian 



2m 



V{x), 



(1) 



where p is the momentum operator —ihd/dx, g is a, non- 
Hermitian field constant in time and space, and V{x) is 
a random potential. Its lattice version is given by the 
matrix 



i) + Vx5x,x'^ (2) 



where x and x' here are site indices, T4; is a random 
potential, and g — ga/Ti with a denoting the lattice spac- 
ing. For simplicity, we focus on the one-dimensional case 
throughout this paper. Periodic boundary conditions 
arc imposed except where stated otherwise. The above 
Hamiltonian reduces to the Anderson localization prob- 
lem for 5 = 0; in this case, it is widely believed that all 
eigenfunctions are localized in one and two dimensions. 

We showed j^,^ that eigenvalues become complex pair 
by pair once g is increased beyond a threshold value 
g = 5ci and that the states with complex eigenvalues are 
delocalized. To show the delocalization, we presented 
two pieces of evidence. First, we numerically showed 
that the states with complex eigenvalues carry a current. 
The current carried by the nth eigenstate is defined by 
jn — den{g) / d{ig), where £„ is the eigenvalue. This is the 
standard definition of the current, because g in Eqs. (0) 
and (^ plays a role of imaginary vector potential. The 
current was clearly nonzero for states in the bubble of 
complex eigenvalues in the band center (see Fig. 13 of 
Ref. and Fig. ^ (b) below), indicating the delocaliza- 
tion of the states. 

As a second indication of delocalization, we showed 
that the delocalized states have complex eigenvalues for 
systems with periodic boundaries, but that all eigenval- 
ues remain real when the same system has open bound- 
ary conditions. This sensitivity to boundary conditions 
is another indication that the corresponding wave func- 
tions are delocalized. We confirmed these two signatures 
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of delocalization in a sufficiently strong imaginary vector 
potential with numerical work and analytic calculations 
on localized impurities. 

This delocalization phenomenon is equivalent to flux- 
line depinning in type-II superconductors with extended 
defects. Suppose that a superconductor has columnar de- 
fects randomly located but mutually parallel and that an 
external magnetic field forces a flux line into the super- 
conductor. The flux line tends to be pinned by a colum- 
nar defect (or a collection of them) when the external 
field is parallel to the defects. When the field is tilted 
away from the axis of the defects, we expect flux-line de- 
pinning at a certain tilt angle (Fig. [l]); see Refs. p^ , p3| , p4[ 
for experiments. 




FIG. 1. Vortex- line system characterized by a one- 
dimensional periodic non-Hermitian transfer matrix. A 
magnetic field H\\ forces a flux line into a cylindrical shell 
of type-II superconductor with columnar defects. The 
current threading the cylinder generates the magnetic- 
field component -ffj^, which tries to tilt the fiux line. 
When this flux-line system is mapped onto a ring of the 
non-Hermitian system, Hx^ becomes proportional to the 
non-Hermitian field g. Below a certain strength of H±^ 
(or g), the fiux line is pinned by a columnar defect and 
forced to run parallel to the defects (the transverse Meiss- 
ner effect [^), except for its slight defiection from the 
pinning center near the top and the bottom of the cylin- 
der For large enough H^, however, the fiux line is 
depinned and wraps around the cylinder as is shown here 
as a helix; This gives rise to a nonzero current that circu- 
lates around the ring of the corresponding non-Hermitian 
system. 



The physics of vortex matter can be mapped onto 
quantum systems with one less dimension by the inverse 
of the Feynman path integral mapping p^]; that is, we 
regard the Boltzmann weight of the fiux line as an expo- 
nentiated action of the world line of a quantum particle 



and make the identification h < — > T, where T is the 
temperature of the vortex system. This procedure gives 
Hamiltonians of the above type. The component of the 
external magnetic field perpendicular to the columns is 
proportional to the non-Hermitian field g. Depinning of 
the fiux line by tilting the field beyond a certain strength 
of g leads to a nonzero current in the corresponding quan- 
tum state 

III. NON-HERMITIAN DELOCALIZATION: 
WAVE FUNCTIONS 

Delocalization of the eigenfunctions themselves was 
not studied directly in Refs. [0J^. The main purpose 
of the present paper is to address this issue. 

Silvestrov computed the right-eigenfunctions associ- 
ated with the model (||) for a 300-site lattice in the re- 
gion of complex eigenvalues and found that their squared 
modulii have a sparse set of well-separated peaks, quite 
different from a conventional delocalized state . How- 
ever, as shown in Ref. ||^, and acknowledged later by 
Silvestrov |2^, it is the product of left- and right- 
eigenvectors which determines the probability distribu- 
tion for a tilted vortex line interacting with columnar de- 
fects deep within the sample. It is this product which 
clearly delocalizes in the conventional sense when the 
eigenvalues become complex. In the hope of avoiding 
further confusion, we first summarize in this section the 
basic relation between left- and right- eigenvectors. We 
then illustrate the delocalization of their product with 
numerical examples from our own extensive 1000-site- 
lattice computations. 



A. Left- and right-eigenfunctions 

We work for concreteness with the continuum Hamil- 
tonian (|l|), but the results also apply to lattice non- 
Hermitian models like (j^). Suppose we have computed a 
set {4)^{x\gy\ of right-eigenfunctions of 'H{g) which sat- 
isfy 

W(g) = e„(g) |n;5) , (3) 

where we adopt the Dirac bra-ket notation, 

^^{x;g)~-^\n-g). (4) 

Although left-eigenvectors need not be simply related to 
right-eigenvectors in general, there is a particularly sim- 
ple relation for the Hamiltonian ([^), which arises due to 
the symmetry [|] 

n\g)^n{~g). (5) 

where f denotes the usual Hermitian conjugate. Indeed, 
as shown below, left-eigenvectors can be obtained from 
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right-eigenvectors by complex conjugation and letting 

9 — ' -g, 

(6) 

or in Dirac notation, 

<t)n{x;g) — > = I"-; (7) 

where f again denotes conventional Hermitian conjuga- 
tion. Our convention that the left-eigenvector {n]g\ is 
defined to be the Hermitian conjugate of |n; — g), not of 
|n;(7), allows manipulations which parallel closely those 
of conventional quantum mechanics. To see that (71; g\ is 
in fact a left-eigenfunction, we calculate 

{n-g\n{g)^{'H{g)^\n--g))^ 

= {n{-g)\n--g))^ 

= (n;5|£„(-5)*- (8) 

Evidently, {n]g\ will indeed be a left-eigenfunction with 
the same eigenvalue as \n]g), provided 

Sni-g)* = £n(5) (9) 

To prove Eq. (||), we let 'H{g) act to the right and the 
left in the matrix element 

{m-g\n{g)\n-g) (10) 

and obtain 

[£m(-ff)* - £n{g)] {m; g\n;g) ^0. (11) 

Equation (^ follows by setting m = n, provided 
{n;g\n;g) ^ 0. More generally, Eq. ( [Tl] ) can be used 
to show that, with proper normalization, the right- and 
left-eigenfunctions form a biorthogonal set, 

(m;5|n;5) = (5„,„. (12) 

This set has the usual completeness relation 

Y^\n-g){n-g\ = \. (13) 

n 

Equation reduces to the usual Hermitian constraint 
of real eigenvalues when (7 = 0. 

Once the eigenvectors are properly normalized, the 
imaginary-time particle propagator is given by 

G(r)=^|n;g)(n;g|e-^'•(^)-/^ (14) 

or, in the coordinate representation, 

G(x,x';r) = ^</,,^(:r)0^(x')e-^"(«)^/^ (15) 

n 

The density distribution of a particle in the ground state 
(which dominates as r — > cxd) is hence the product of left- 
and right-eigenfunctions, 0gs(a;)0|^(x). As was shown 
in Ref. this product gives the probability distribu- 
tion of a flux line far from the sample boundaries in the 
imaginary-time direction (the top and the bottom edges 
of the cyhnder in Fig. |^). The square moduhi |0^(a;)p 
and |0^(a;)p are irrelevant for the bulk properties. 



B. Delocalization of <\>^<\>' 



We now illustrate the different behaviors of |(^^(x)p, 
|(/)^(a:)p and (j)^{x)(j)'^{x) with numerical results, empha- 
sizing that the product ^^</>^ is clearly delocalized in 
the conventional sense when eigenvalues become com- 
plex. We consider a particular realization of the random 
Hamiltonian (||) on a 1000-site lattice. The parameters 
are set to t = 2 and g = OA with the value of Vx at each 
site chosen randomly from the range [—1.5, 1.5]. These 
values are the same as used in Refs. except that 

the system size is greater in our calculation. (Note that 
the definition of t differs by factor two.) The energy spec- 
trum is shown in Fig. ^ (a). The states between the two 
mobility edges Ec — ±2.34 have complex eigenvalues (and 
hence we would argue are delocalized), while the other 
states are localized. Every delocalized states carries a 
complex current as is shown in Fig. |^ (b). The imagi- 
nary part of the current determines the tilt angle of a 
flux line 

Figure^ (a) shows the functions (j>^^{x)<j>^{x), \<j>^{x)\'^ 
and |0^(a;)p for the (localized) ground state. All quanti- 
ties are normalized so that the summation over x yields 
unity. We stress, however, that the normalization only 
makes physical sense for (j)^[x)<f)^{x). Everywhere in the 
regime of localized states, 

0«(x;g)cxe^^/''(/.„(a;;O) 
and (/.^(x; g) cx e"^^/" (0„(x; 0))* (16) 



for large enough systems, where <j)n{x]Q) is the wave 
function of the Hamiltonian with g — 0. Note that 
Eq. (H) is obeyed. (The specific g-dependence in Eq. (|l^ ) 
only holds for localized states, for which we can al- 
ways choose 0„(x;O) to be real.) Hence, the product 
<j)^^{x; g)(j)^{x] g) = |(/)„(a;;0)p does not depend on g un- 
til the state is delocalized for large enough g. This is a 
mathematical expression of the transverse Meissner ef- 
fect, or the rigidity of the pinned flux line against the tilt 
of the applied magnetic field; see Refs. P,p5[ for details. 

Figure |3| (b) shows the functions |0^(a;)(/)^(a::)|, 
\(f)'^{x)\^ and |<^^(a;)p for a state slightly below the lower 
mobility edge, while Fig. || (c) shows those for a state 
slightly above the edge. (We only plot the amplitude of 
the function (j)^{x)(t)!^{x); the phase oscillates rapidly for 
delocalized states away from the band edges.) The func- 
tion \(j)f^{x)(j)^'{x)\ changes dramatically across the mobil- 
ity edge, while the changes in |0^(a;)p and |0^(a;)p are 
less noticeable. 

The delocalized nature of (j)^{x)4>^{x) appears even 
more dramatically deep inside the bubble of complex 
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eigenvalues. Figure |^ (a) shows the same functions 
for an eigenstate with the eigenvalue e„ = —2.01239 + 
i 0.200376. This is the 166th state, which roughly cor- 
responds to the 50th state of the 300-site system stud- 
ied by Silvestrov ^j^. In Fig. | (a), the function 
|(^^(a;)0^(x)| is extended and approximately constant, 
while and exhibit a sparse set of well- 

separated maxima. Following Silvestrov, we plot the log- 
arithm of these functions in Fig. ^ (b). The product of 
(f>^{x) and (j)^{x) is remarkably constant and is extended 
in conventional sense. On the other hand, the ragged 
wandering nature of In and ln|0^(a;)| is consis- 

tent with the conjecture pO| , pT[ that these functions be- 
have like random walks as a function of x; this is the 
subject of the next section. 
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FIG. 2. The complex energy spectrum and the current 
distribution of the Hamiltonian (^) on a 1000-site lat- 
tice. Each eigenvalue is marked by a tiny cross in (a). 
Each pair of complex conjugate eigenvalues in (a) has 
the current shown in (b), with the reals part of the oppo- 
site sign (tiny crosses) and the identical imaginary parts 
(the dashed line). The parameters in Eq. are set to 
t — 2 and g = 0.4 with each Vx chosen randomly from 
the range [—1.5, 1.5]. 
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FIG. 3. The functions |(/)^(x)0-"(x)| (thick solid lines), 
\4)^{x)\'^ (dashed lines) and |(/)"^(a;)p (dotted lines) for 
the following cases: (a) the ground state (e = —2.94682); 
(b) an eigenstate just below the lower mobility edge 
(the 72nd state with e = —2.35214); (c) an eigenstate 
just above the lower mobility edge (the 80th state with 
£ = -2.31828 -t- i 0.0103479). The system is the same 
as the one used in Fig. ||. The serial number of each 
state represents the ascending order of the real part of 
the eigenvalue. 
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FIG. 4. The functions (thick solid hnes), 

|0^(a:)|^ (dashed lines) and |0^(x)p (dotted hnes) for the 
an eigenstate in the delocalized regime (the 166th state 
with e = -2.01239 + ^0.200376); (a) a linear plot and (b) 
a semi-logarithmic plot. The system is the same as the 
one used in Fig. 0. 



IV. RANDOM- WALK BEHAVIOR OF 
EIGENFUNCTIONS 



In this section, we turn our attention to the left- and 
right-eigenfunctions considered separately for large g. 
One of the interesting results of Silvestrov [^,0 is a 
random-walk-like behavior hidden in the logarithms of 
the modulii of these eigenfunctions. We first illustrate the 
random-walk behavior with our more extensive numeri- 
cal results and then show that, at least for the ground 
state, Silvestrov's observation is related to earlier results 
obtained for charge density waves |^ and population bi- 
ology ||l2|. As a concluding remark, we argue that the 
behavior of sample-to-sample fluctuations of the left and 
right ground-state eigenfunctions considered separately 
for large g is "intermediate" between that expected for 
localized and delocalized states. 



A. Vortex-line distribution at boundaries 



What information is contained in the functions 4'n{x) 
and (pni^) (considered separately) for flux line systems? 
For a single vortex line, only the ground state contributes 
in the limit of a very long cylinder. As discussed in 
Ref. PI a nd exploited in a very recent paper by Silve- 
strov |26| , the (nodeless) ground state wave functions 
<j>gs{^) and (j)gs{x) {not their modulii squared) are pro- 
portional to the vortex-line probability distribution at 
the boundaries where it enters and leaves the cylinder 
(see Fig. I). 

There are then two cases to consider. For small and in- 
termediate values of 5, the spectrum is either completely 
localized or only partially delocalized as in Fig. § (a). In 
this case the ground state is localized as in Fig. || (a), 
and hence the single vortex line is pinned close to a pre- 
ferred columnar defect in the bulk of the superconductor 
cylinder. The right- and left-eigenfunctions are shifted 
relative to their product. This reflects the tendency of 
the localized vortex line to tear away from the pinning 
center at the top and bottom of the sample when g is 
nonzero; see Fig. 15 (a) of Ref. |^ for a demonstration. 

The second more interesting case is for large g > 
gc2, such that all states, including the ground state, 
are delocalized. Using the WKB approximation, Sil- 
vestrov 1^ , ^ argued for random-walk behavior of the 
logarithm of the wave functions in this case. For con- 
creteness, we show some of our numerical results for the 
ground state of a 2000-site lattice. (Silvestrov |2^ did 
not show numerical results for the ground state in this 
regime.) Figure || (a) shows the ground state quantities 
(/)^s(x), and (j)^^{x)(j)^^{x) for g = 1.5h/a > 5^2; the 

values of the other parameters are the same as in the ear- 
lier figures. The product (/)gg(x)(/)^(x) is approximately 
constant, while the (nodeless) eigenfunctions (j)gsi^) ^^'^ 
4>gsi^) are quite different than in Fig. || (a): They exhibit 
multiple sharp maxima which are rather well separated. 
In view of these multiple maxima, one might question 
whether it is appropriate to call such eigenfunctions "lo- 
calized" [^0. Figure || (b) shows the same ground 
state quantities in a semi- logarithmic plot. The wander- 
ing, ragged shape of ln0j(x) and \n(j)gg{x) indicates the 
random- walk behavior. 

The different shapes of (?!)gg(x)(/)J(a;), 0j(a;) and 0gs(a;) 
reflect the different optimization problems of the vortex- 
linc configuration in the bulk, at the top and at the bot- 
tom of the superconductor cylinder (see Fig. Since 
the string tension of the vortex line (the "mass" of the 
corresponding quantum particle) is missing outside the 
superconductor, the vortex line can take better advan- 
tage of the potential energy at the top and bottom of 
the sample than in the bulk; hence the sharp maxima in 
(/)gg(a;) and (j)gg{x). The multiple maxima indicate that 
the depinned vortex line can enter and exit the super- 
conductor at variety of preferred locations. 
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The optimization problems are also different at the top 
and the bottom of the cylinder in Fig. |l]. When the vor- 
tex line enters the sample from below, it is the succession 
of defects counterclockwise to the entry point which are 
most important. When exiting the sample, it is the de- 
fects clockwise to the exit point that matter most. Hence 
the peaks in (f>gs{x) and in 0j(a;) appear at very differ- 
ent locations. Nevertheless, the entry and exit probabil- 
ity distributions are strongly correlated with each other, 
since their product is approximately constant. 



B. Renormalization group for the ground state 

In the following, let us reproduce Silvestrov's WKB 
result for delocalized states in a more controlled approxi- 
mation. For the ground state, his result is in fact a special 
case of numerical, scaling and renormalization group cal- 
culations applied previously to related problems in one- 
dimensional charge density waves |^ and population bi- 
ology in d dimensions . 

We start with the time-dependent Schrodinger equa- 
tion for the continuum Hamiltonian (PO), 
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-V{x)^Ij"{x,t). 



(17) 



We assume uncorrelated finite-width randomness of the 
potential. 



V{x)V{x') ^/^^5{x-x'), 



(18) 



where the overbar denotes the random average and A is 
the width of the random distribution. 

The d-dimensional generalization of Eq. ( |l7|) was stud- 
ied in Refs. ||^,|l^ via the "Cole-Hopf transformation," 



tp'^'{x, t) — exp 



^{x,t) 



2mh 



(19) 



The "Cole-Hopf transformation" is just another name for 
the WKB method. The second term in the exponent of 
Eq. (|l|) is added in order to offset the ground-state en- 
ergy (which has no effect in the physics of flux line) . The 
equation for <i> is 

dr m dx 2m dx"^ 2m \ dx J 

To see the relevance of each term in the long-distance 
limit, we change the scale as part of a renormalization- 
group calculation, according to 



(21) 
(22) 
(23) 



FIG. 5. The functions (p^ {x)(f>"{x) (think sohd hues), 
(j)^{x) (dashed lines) and ^f)^(x) (dotted lines) for the de- 
localized ground state of a 2000-site lattice; (a) a linear 
plot and (b) a semi-logarithmic plot. All these functions 
are positive definite in this case. The normalization of 
(j)^{x) and 4>^{x) is different from the one in earlier fig- 
ures; each is normalized so that its sum over x (not the 
sum of the squared modulus) becomes unity. The param- 
eter values are the same as the ones used in Fig. |^ except 
that g — 1.5. 



where the exponents z and a are determined below. Thus 
we have 



5$ 



m dx 
2m dx^ 



z+a-2 



1 



2m \ dx 



(24) 



The rescaled random potential is defined by V{x) = 
si/2y(sj) so that it satisfies V{x)V{x') = A'^S{S: - x'). 
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The first term of the right-hand side of Eq. (g4|) is a 
drift term and the second term is the random potential 
term. To keep these two terms fixed in the long-distance 
limit s — *■ oo, we set z — 1 and a = 1/2. The third 
and the fourth terms are then irrelevant variables in a 
perturbative renormalization group like that constructed 
in Ref. |12|. Thus a Gaussian fixed point controls the 
physics of what later turns out to be the regime g > gc2- 

Upon defining renormalized parameters by 



TO = S ^ "+^TO, 



and 



(25) 
(26) 
(27) 
(28) 



we arrive at a Langcvin-type equation in the long- 
distance limit, 



OT TO OX ' 



(29) 



where we have dropped the "tilde" symbol from all quan- 
tities. Since the ground-state energy was already offset in 
Eq. (19), we can eliminate the time derivative by moving 
onto a new set of coordinate as (x,r) — > {x,T~{m/g)x). 
We thus see that the solution is a random walk evolving 
into the x direction: 



(a) 



(b) 



(c) 
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^(x,t) = ^ix) = — / V{x')dx 
9 



(30) 



The stationary right-eigenfunction in the long-distance 
limit is hence given by 



4>^,{x) = exp 



.HL I vix')dx' 

gh 



(31) 



except a normalization factor. This is equivalent to the 
ground-state (fc = 0) solution of Silvestrov's calcula- 
tions |2|,|2^,|§. Equations (§) and (H) then give the 
left-eigenfunction as 



exp 



^ r v{x')dx' 

.9^ J 



(32) 



Note that the random-walk behavior disappears for the 
product (/)gg0j, which is completely homogeneous in this 
approximation. This is consistent with the numerical re- 
sult in Fig. I for g > gc2- 



C. Sample-to-sample fluctuations 



Silvestrov [ poy21|j26| referred to the above behavior of 
(f)^^ and (/)gs as "stochastic localization" (or simply as "lo- 
calization" in some sentences). Although the behavior is 
quite different than the smooth delocalized behavior of 
their product, it is not clear to us whether such states 
should be called "localized" either. 



FIG. 6. Schematic views of sample-to-sample fluctu- 
ation W{Lx) for various types of ground state: (a) 
"random- walk" wave functions of the form (31); (b) a 



set of wave functions localized in a conventional sense; 
(c) wave functions which are extended in a conventional 
sense. Different curves in each graph indicate wave func- 
tions of different samples. 



To stress this point further, we follow Refs. |,|l| and 
consider the sample-to-sample fluctuations of the loga- 
rithm of the ground-state wave function at a flxed loca- 
tion. 



W{L^) = \ii(f)gsix)^ - ln?!)gs(a;) 



(33) 



where is the system size in the x-direction and 
(f>gs{x) is either the right- or the left-eigenfunction of 
the ground state. (We choose a normalization such that 
/ 0gg(/)^(ia; ~ 1.) The a:-dependence of this quantity 
should disappear owing to the statistical translational in- 
variance. As is shown below and in Fig. |6[ we would have 

• W(L^) ~ 0(L.^) for the ground-state wave func- 
tions d^) and (11); 

• W{Lx) — 0{L^) for a conventional (Hermitian) lo- 
calized ground state; 



ground state. 



for a conventional extended 
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From this point of view, the random- walk behavior of 0^ 
and may be viewed "intermediate" between locahzed 
and delocalized. 

The first La;-dependence of the quantity WjLx) is de- 
rived either from the wave function (pi]) or ([3^) 



as 



W{Lx) oc // dx'dx"V{x')V{x") = 0[Lx). 



(34) 



Next, to calculate W{Lx) for a conventional localized 
ground state, we assume its asymptotic form as (j){x) ~ 
eicp{—K\x — Xc\) and that the value of k is approximately 
equal for all samples but the localization center Xc is dif- 
ferent from sample to sample; see Fig. |^ (b). The random 
average in Eq. (|3|) then reduces to the average over Xc- 
Thus we have 



WiLx) 



Lx 
0{Ll). 



dXr 



K 


f 


Tx) 


1 X X Q 1 dx Q 



(35) 



Finally, the logarithm of a conventional extended ground 
state should be approximately homogeneous in space for 
all samples and hence have little sample-to-sample fiuc- 
tuation as illustrated in Fig. |^ (c). This is the behavior 
of when g > gc2- 



V. SUMMARY 

In conclusion, we have argued that, contrary to some 
recent statements in the literature, delocalization does 
appear in the eigenfunctions of Hamiltonians such as (|^) 
and (^ when the spectrum becomes complex, provided 
one studies the product of left- and right-eigenfunctions. 
Delocalization defined by this criterion is consistent with 
earlier definitions based on ability of states with complex 
eigenvalues to carry a nonzero current and the sensitivity 
to boundary conditions 0J^. 

The left- and right-eigenfunctions considered sepa- 
rately provide interesting information about the physics 
of an isolated vortex line at its entry and exit points. 
Similar conclusions apply to interacting arrays of vor- 
tices. See Sec. VIII of Ref. [|[ for a discussion of this 
non-Hermitian many-body problem. An interesting in- 
vestigation of tilted interacting vortices at the entry and 
exit boundaries has been initiated by Silvestrov [E6|. 
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